    
a = 0;
b = 48;
ya = 10;
h = 0.01;
func = 'f1';
kArr=[0.25,0.5,1,2,5,10,15,20];
global k;

for i=1:length(kArr)
k=kArr(i);
RK=rk4gral(a,b,ya,h,func);
max(RK(:,2))
subvec=RK(round(length(RK(:,1))/4):length(RK(:,1)),2);
amp(i)=max(RK(:,2))-min(subvec);
plot(k, amp(i), 'r*');
hold on;
end

plot(kArr, amp, 'r-');

legend('Amplitud de la Temperatura', 'Location', 'NorthEastOutside');

%title('Amplitud de la temperatura para distintos valores de k')
xlabel('k');
ylabel('Temperatura (C)');

hold off